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Abstract. - We present results for a model that describes a quantum point contact. We 
show how electron-electron correlations, within the unrestricted Hartree-Fock approximation, 
generate a magnetic moment in the point contact. Having characterized the magnetic struc- 
ture of the contact, we map the problem onto a simple one-channel model and calculate the 
temperature dependence of the conductance for different gate voltages. Our results are in good 
agreement with experimental results obtained in GaAs devices and support the idea of Kondo 
effect in these systems. 



Ballistic transport in quantum wires with ideal contacts leads naturally to conductance 
quantization in units of 2e'^/h [1-3]. This notable result has become one of the paradigms 
of mesoscopic physics. By a quantum wire we understand conductor that is so thin that the 
quantization of the transverse modes becomes important. In short quantum wires we may 
also expect some quantization effects along the wire axis [4] . In some of these short wires and 
in quantum point contacts, as the gate voltage is varied, together with the 2e^//i plateaus the 
conductance shows some structure near 0.7(2e^//i) [5-8]. This structure has been systemati- 
cally observed GaAs devices and is known as the 0.7 anomaly. Recent experiments presented 
evidence on the magnetic nature of this structure: for gate voltages at the 0.7 anomaly the 
conductance increases as the temperature is lowered, an effect that may be interpreted as 
a signature of Kondo effect [7]. The Kondo picture is supported by the scaling of the ex- 
perimental data, that close to the 0.7 anomaly shows a universal temperature dependence. 
Moreover the low temperature differential conductance presents a narrow structure at zero 
bias voltage, which splits with an in-plane magnetic field as due to a Kondo resonance. The 
conventional Kondo effect is the magnetic screening of a localized moment [9], if this effect 
is taking place in the quantum point contact (QPC), then the main question is: where and 
how is the magnetic moment generated? Recent spin dependent functional density calcula- 
tions suggest the formation of a magnetic moment localized at the QPC for a device built 
on GaAs [10-12]. This scenario in which a magnetic moment is spontaneously generated by 
the electron-electron correlations is appealing since it naturally leads to the Kondo physics 
observed experimentally, however there are still many open questions related to the magnetic 
nature of the QPC and its relation to the 0.7 anomaly. 
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Fig. 1 - (a) Schematic picture of tiie model for the QPC, the potential at the constriction is A. (b) 
Magnetic moment (continuous line) and total charge (dashed line) at the constriction of a 3-site- 
wide and 9-site-long QPC. The first conduction channel opens at a and the second at b, the large 
magnetization at these points signals the formation of a local magnetic moment. The local interaction 
isU = 2.0 and the Fermi energy is ep ~ —0.95, with all parameters in units of t = 1. The source and 
drain slabs have a width of 10 sites, (c) Detail of the total charge around a (dashed line) compared 
with a non- interacting system (continuous line). 



For the sake of simplicity we resort to a one band Hubbard model with local interactions. 

i,(y i {i,j),a- 

here Ci^ annihilates an electron with spin a at site i, U is the on-site Coulomb repulsion and 
the last term in Eq. describes the nearest neighbor hopping with matrix element t. The 
on-site energy Si describes the potential at the QPC. For simplicity, in what follows we use 
a square potential parametrized with a single variable A: using the geometry of Fig. Ql^a) 
we take = for i in the source and drain, and Ei = A, oo for i in the contact and at 
the sides of the contact respectively. The energy A can be controlled by a gate potential 
Vg and we may consider A oc — |e|V^. The square potential maximizes resonances for states 
with wavelengths commensurate with the point contact, nevertheless smoother potentials 
give similar results [13]. We look for the ground state of Hamiltonian (1) in the unrestricted 
Hartree-Fock approximation for slabs with a QPC as in Fig. da). In the central region 
[between dashed lines in Fig. ^a)], the expectation value of the electron number for each 
spin and at each site is calculated self-consistently. Away from the constriction, in the source 
and drain, we use the Hartree-Fock solution of a uniform system with approximately 0.5 
electrons per site. For the results presented below, we checked that the boundary between 
the two regions is far enough from the contact so that the solution becomes independent 
of its position. Given the geometry, i.e., the length and the width of the point contact, 
and the value of U within a wide interval, we find that for some values of A a magnetic 
moment is formed in the region of the constriction. We define the total charge Q and the 
magnetization A4 of the contact as the expectation value of the operators Q = ^ cJ^Cicr 

and M — ^ J2i ~ '^I'^ilj where the sum is over all the sites in the neck of the QPC 

[shaded region in Fig. Qla)]. The results are summarized in Fig.QJb): for large negative gate 
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voltages the total charge Q in the point contact is exponentially small. As A decreases Q 
increases and for some values of A there is an abrupt increase in the total charge. At the 
step-like structures obtained around the point marked as a in Fig.^b) [see also detail in^c)], 
approximately one electron is transferred from the source and drain to the point contact, a 
behavior characteristic of Coulomb blockade. Between the two first steps in Q there is a 
spin 1/2 localized at the point contact as indicated by the magnetization curve shown also 
in Fig. ^b). In Fig. ^c) we present a detail of the behavior of the charge in the QPC as a 
function of the gate voltage near the point a of Fig. Q^b) for two values of the local interaction. 
For the U = 2 case the entrance of an electron in the QPC is blocked, between the two steps, 
by the Coulomb interaction. If U is decreased the separation between the steps diminishes 
linearly, and also the magnetization decreases as the number of electrons in the QPC is less 
well defined. For U < 0.85 the two steps merge into a SQ ^ 2 step and the magnetization 
vanishes. 

At the point a of Fig. ^b), where a spin 1/2 local moment is stable, a narrow resonance 
is crossing the Fermi energy. This resonance is clearly seen in the local density of states 
(LDOS) of the QPC as shown in Fig.|2Ia). This first resonance corresponds to a state with 
no nodes in the QPC as can be inferred from the magnetization of Fig. |2Ic). Note that 
the magnetic moment is localized at the QPC. At the point marked as b in Fig. ^h) again 
a spin 1/2 is localized in the QPC, this is due to a resonance of the second channel [see 
Fig.|2Ib)] that generates a magnetization profile with a node along the mayor axis of the QPC 
as shown in Fig. |2Id). When the chemical potential coincides with a narrow resonance, the 
system Hamiltonian can be approximately mapped into that of an Anderson impurity [14] 
with effective interaction Ueff and resonance width Fe//. If Ueff ^ TrFe//, in analogy with 
Anderson criterion for magnetic impurities [9], we expect a magnetic moment in the QPC. 
The effective parameters scale with the size of the QPC as: Ueff ^ U/A where A is the area 
of the QPC and, for the geometry of Fig. QJa), Fg// oc where L is the length of the 

QPC [13]. 

These results deserve a few comments: 

i.- The behavior of Q as function of A and the occurrence of a local magnetic moment in the 
QPC is reproduced for a wide range of values of the Coulomb repulsion U. For the electron 
densities used up to now (~ 0.5 electrons per site), the critical value of U that generates 
a global instability Uc is larger than 8t. The behavior presented in Fig. ^b) corresponds 
to U = 2t, far away from the global instability. We stress that the model aims to give a 
qualitative description with the minimum number of parameters. 

a. -The formation of a local magnetic moment is not a peculiarity of the model used to 
simulate the potential due to the gate voltage. We have analyzed different potentials to 
show that the local instabilities shown in Fig. ^b) , are robust for potentials that have a flat 
region and produce resonances. Conversely, for wedge-like point contacts [1], the Hartree-Fock 
solution does not show a magnetic moment for moderate values of U. 

Hi.- For the geometries of Fig.^a), narrow resonances are obtained at the bottom of each 
channel band. Localized magnetic moments are then obtained each time a new channel be- 
comes active in agreement with Ref. [8] . The details of the magnitude and spatial distribution 
of the magnetic moment depend on the shape of the QPC potential. 

iv.- The Hartree-Fock approximation, as well as other methods like the spin dependent 
functional density approximation [12], can not be fully satisfactory since they break the lo- 
cal symmetry. In the exact solution of the problem, spin fluctuations should recover local 
rotational invariance. However as for the original problem of magnetic moment formation 
in impurities [9], the Hartree-Fock results do indicate the region in parameter space where 
we can expect magnetic fluctuations to play a central role in the low energy physics. A de- 
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Fig. 2 - (a) Average in the QPC region of the total local density of states (spin up plus spin down) for 
the system of Fig.IHb) and A = 2.27. (b) Same as (a) with A — 0.66. The vertical dashed line in (a) 
and (b) marks the position of the Fermi level, (c) and (d) spatial distribution of the magnetization 
for the (a) and (b) situations respectively. 



tailed analysis of the local magnetic instability for different QPC potentials will be presented 
elsewhere [13]. 

To go beyond the Hartree-Fock approximation and to make a connection with the Kondo 
problem a simplification of the model is required. To do so, we first revisit the transport 
properties of a non-interacting QPC that can be calculated exactly using standard techniques 
[3,15]. The conductance as a function of the gate voltage Vg is shown in Fig.Ofa). The results 
presented in the figure are characteristic of potentials having a flat region along the mayor 
QPC axis, the details however depend on the specific potential used to describe de QPC. The 
conductance for longer QPC oscillates more rapidly since resonances are narrower and closer 
in energy. We note that here, when the conductance in units of 2e^//i is smaller than one, 
the transmission is due to a single channel. In this one channel regime, and at low enough 
temperatures, we can map the problem onto an effective one-dimensional (ID) model. 

In what follows we present a model that, in the gate voltage interval indicated with vertical 
dashed lines in Fig. Eta), reproduces the structure of the conductance and of the density 
of states. In Fig. Elb) a schematic picture of a linear chain with a 3-site central region 
representing the QPC is shown. The Hamiltonian is given by 

oo oo 
— oc.cr i— — oo,(T 

To reproduce the structure we adjust the hybridization of the central site V — to,i — t-i,Oy 
those of the adjacent sites ii_2 — ^-2,-1 and the diagonal energies eo and e_i — ei. Only the 
diagonal energies in the 3-site central region are modified by the gate potential (they are shifted 
by V^). The other paremeters are: ei — for |i| > 1 and U^i+i = 1 for i ^ {—2, —1, 0, 1}. The 
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Fig. 3 - (a) Conductance for a 3 x 9 QPC in a non-interacting system as a function of the barrier height 
in the QPC region, (b) Effective ID model to describe the first step in (a), (c) Conductance and (d) 
local density of states at the central site of (b) for a QPC with eo = ~Vg, ei — e_i — — Vg + 0.15, 
V = 0.15/%/2, and ti,2 = t-i,-2 = 0.25. 



conductance as a function of a gate voltage for this ID model, shown in Fig.|2Ic), reproduces 
the first structures of Fig.lSfa). The first peak in the conductance is due to a resonant state 
crossing the Fermi level. This resonance is observed in the local density of states (LDOS) 
of the central site as shown in Fig. |3Jd). Having this non-interacting ID model that mimics 
the transport properties of the QPC, we can now turn on the interaction that generates a 
local moment. The local moment should be associated with the first resonance in the LDOS 
and then we include a local interaction at the central site: TCid = T^id ~^ U'^oi'^oiCq^cqi. 
The temperature dependence of the conductance for this model in the linear regime can 
be calculated exactly using the numerical renormalization group [16, 17]. We proceed as in 
Ref. [18] to put the conductance in terms of the spectral density of the correlated site po{lu)- 

G{T) = ^27r J dc.(-^)r(c.)po(^) (3) 

where /(w) is the Fermi function, and r(cLi) = TrV'^pi{uj) with pi{oj) the spectral density of 
site 1 with V = 0. 

The results are shown in Fig.^fa) and01b) for two sets of parameters that correspond to a 
long QPC (sharp resonances) and a short QPC (wide resonances) respectively. At intermediate 
temperatures the conductance presents a plateau at about 0.7(2e^//i). This structure is not 
universal, depending on the parameters a more pronounced maximum can be obtained and 
the whole structure may be slightly shifted up or down. As the temperature decreases the 
conductance increases toward 2e^//i, the plateau shifts away from 0.7(2e^//i) and disappears. 
These numerical results are notably similar to the experiments of Ref. [7] and also agree 
with the experiments of Ref. [4] where sharper structures are found for longer QPC and a 
maximum, instead of a plateau, is obtained in the conductance. 

When the first resonance of the QPC is above the Fermi level [Vg < —0.12 for the param- 
eters of Fig. EJa)] the conductance and the charge in the QPC are small. An increase in Vg 
pushes the levels of the QPC below the Fermi level, favoring the charging of the QPC. At 
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Fig. 4 - (a) Conductance as a function of the gate voltage Vg for the one-dimensional model of 
Fig. 2(b) with U = 1.0 in the central site and the other parameters as in Fig. |5Jc). (b) Same as 
(a) with V = Q.2/y/2 and ei = e-i — —Vg + 0.25. (c) Conductance as a function of the scaled 
temperature for the parameters in (a), (d) Kondo temperature obtained from the scaling in (c), 
= TxiVg = -0.025) = 1.35 x 10"^ 



Vg ~ —0.08 the Fermi level of the leads coincides with the first resonance in the density of 
states of the QPC. The charge in the QPC is ~ 1, and it has a spin 1/2. For Vg > —0.07 the 
magnetic moment is formed and at low temperatures there is an increase in the conductance 
as the Kondo effect develops. The conductance as a function of temperature for different 
values of the gate voltage is shown in Fig. 0Jc). At low temperatures, all the curves collapse 
into a single one when the temperature is properly scaled. We define the Kondo temperature 
Tk as the scaling parameter. Contrary to what is obtained in the conventional Anderson 
model [9] , in our mapping the Kondo temperature increases as the gate voltage increases and 
the resonance is pushed away from the Fermi level [see Fig.0Jd)]. This is in agreement with 
the experimental observations and is due to an increase of the density of states at the Fermi 
energy as Vg increases and the second longitudinal resonance approaches the Fermi level. The 
main difference with the conventional Anderson model is the dependence of the hybridization 
function r(w) on Vg. In the Anderson model r(a;) is usually taken as a constant r(ijj) = r(ei?) 
to describe e.g. quantum dots. In this case r(ei;') increases as the gate voltage increases, lead- 
ing to an enhancement of the Kondo temperature. We presented results for U = 1.0, however 
similar results are obtained for U > 0.25. To be able to observe the Kondo behavior the local 
interaction U must be larger than the minimum value of the hybridization r(ei^). 

In summary we have studied the formation of a magnetic moment in quantum point 
contacts and its relation with the observed anomalies in transport experiments. The Hartree- 
Fock picture presented here accounts for the occurrence magnetic moments when resonant 
states of the QPC are aligned with the Fermi level. For a moderate local interaction, the stable 
local magnetic moments are related only to the first longitudinal resonance of each transverse 
mode, i.e., at the beginning of each conductance step, in agreement with the results of Reilly 
et al. [8]. However, the behavior of the conductance at each conductance step will depend of 
the particular details of the QPC potential. These results are robust for QPC that generate a 
flat potential along the contact axis. By mapping the problem onto a one-channel model, we 
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calculated the temperature and gate voltage dependence of the conductance that are in good 
agreement with the experimental observation. Our results support the idea of Kondo effect 
in the point contact however our proposal differs from previous ones: in the work by Meir et 
al. [12] it is assumed that as the gate voltage is increased, the magnetic moment is formed well 
before the conductance increases although it is essentially decoupled from the electron gas. 
In our approach the conductance increase in accompanied by the formation of the magnetic 
moment and this occurs each time a new channel starts to participate in the transport trough 
the point contact. 

A large magnetic field may destroy the Kondo correlations and could be used as an external 
probe to control the spin fluctuations. However, in the experiments of Cronenwett et al. [7] 
the Kondo temperature has been estimated to be as large as 10 K and consequently only 
huge magnetic field could destroy the Kondo correlations for all the gate voltages. In a 
recent publication [19], other GaAs point contacts that for zero magnetic field present the 
0.7 anomaly, for fields of 15 Tesla show a conductance plateau at e^/h with a structure that 
clearly resembles a resonance at the edge of the first conductance step in agreement with the 
prediction of the model presented above. 

* * * 
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